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PREFACE 





The research reported here was done at the U. S. Naval Postyraduate 
School, Monterey, California during the period February through April, 
1954e The work was undertaken to determine experimentally the maximum 
stress and the torsicnal rigidity of a rowd shaft with diametrically 
opposite parallel flat sides. This information was obtained for several 
configurations by direct experiment and for one configuration by an 
electrical analogy. 

The author is greatly indebted to Prof. R. E. Newton for many 
helpful suggestions and much constructive criticism and assistance in 
the preparation of this manuscript. In addition the author wishes to 
thank Mr. Norman ‘ialker for the extreme care and patience exercised 
in the preparation of the test specimens. It is further desired to 
acknowledge the help of F. N. Lauck, Ph3, USN, who did the photographic 


work which appears in the thesis. 
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This thesis describes and reports the results of an experimental 


determination of the maximum shearing stress and twisting moment of cir- 
cular eee with diametrically opposite flat sides. Five AL-61 Aluminum 
alloy specimens were tested, one fully round shaft and the remaining four 
with progressively larger flats. Strains were measured with resistance 
wire strain gages. Angles of twist were determined by optical means. 

An electrical analog experiment was conducted to determine independ- 
ent results for one of the configurations. 

The experimental results obtained by both techniques compared very 


favorably with mathematical solutions obtained by previous investigators. 














CHAPTER I 


AN INTRODUCTIGN TO THE PROBLEM 


1. Statement of the problem. 

This thesis is a report of an experimental determination of the 
torsional rigidity and maximum torsional stress occurring in a solid 
round bar having diametrically opposite flat sides. 

2. State of stress for a solid round shaft with diametrically opposite 
parallel flats. 

The state of stress can best be visualized by reference to the 
membrane analogy for the torsion problem {1 7. Considering the de- 
flected surface of the membrane conclusions can be made as to the stress 
distribution. The shear stress in any direction is proportional to the 
slope of the surface perpendicular to that direction. It varies from 
a@ maximum tangent to the contours to zero perpendicular to the contours. 
Thus the maximum shear stress acts at points where the contour lines are 
closest to each other. With this in view reference to Fig. 1, which is 
@ cross section of the shaft, leads 
one to conclude that the maximum 
stress will be at the center of 
the flat sides and zero at the 
center of the cross section and 
at the four corners. 

Fig. 1. Cross section of 
the shaft. 
3. Previous solutions to the problem. 


H. Okubo [ 2_/obtained an approximate solution for maximum shearing 











stress employing a variation of the technique given by Shepherd / 3G 
Okubo obtained an approximate solution in the form of elementary functions 
which satisfied the differential equation but only approximately satisfies 
the boundary conditions. W. J. Carter and J. B. Oliphint / 4_/ obtained 
solutions for torsional rigidity and maximum shear stress by numerical 
relaxation. For very small flats Carter and Oliphint made use of a 
hydrodynamic analogy to obtain the maximum stress. The above authors 
carried out the computations for several configurations, but no experi- 
mental results were presented. A search of the available literature by 
the writer disclosed no further work on the problem. 
4. Electrical Analogy 

The problem considered is that of twisting a prismatical bar with 
couples applied at the end cross sections (St.-Venant torsion Lome 


It is assumed that the bar is homogeneous, isotropic and obeys Hooke's 


law. y 





Fig. 2. Bar in Torsion. 


Ste-Venant showed that the only stress components different from 
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Fig. 3- Cross Section of the Shaft. 


The stress function satisfies the partial differential equation 
2 . : 
V*o= $e + $4 =-26¢ A ae 


and the boundary caomdition 


= O on C. (33 


G is the shear modulus for the material and (> is the anzle of twist per 


unit length. R is the region of the cross-section of the shaft and C is 


= 




















its boundary curve (see Fig. 3.) The twisting moment is found by 


Mall g ax dy (4) 
PR 


To transform the problem, assume the stress function has the form 
g= 9 -— SP Oty) (5) 


in the region R of Fig. 3 and is equal to zero on the boundary Ce. Put- 


ting this in Eq. (2) gives 


Vo =o (6) 


which is Laplace's equation to be satisfied in Re. The bowdary condition 


now reduces to 





p= BE (42 +y?) on C. (1) 


Thus the torsional problem reduces to obtaining a solution of Laplace's 
equation over the cross section and satisfying the boundary condition 
Eq. (7). 

There are many techniques available for obtaining a solution of 
Laplace's equation. The author chose to use a General Electric Analog 
Field Plotter with Teledeltos uniform conducting paper because of its 
simplicity, ease of use and the availability of components. The experi- 


mental assembly used is shown in Fige 4- 
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gp =K(V) 


where K is a constant of proportionality and V is the potential. Then 
by Eq. (7) the potential varies parabolically along the flat side and is 
constant along the arc. The equipotential lines are lines of comstant 
p’. From the equipotential lines and K, the Q’ contours can be obtained 
for the region. Knowing on enables one to obtain the gD contours by 
means of Eq. (5). 

Eqse (1) now become, by substitution of Eqs. (5) and (8), 


Tig=htt -Gey ,  Typ=-k¥ + GOx 7 


where 7 and . are the components of the potential gradient as 


obtained from the equipotential lines. Eq. (4) now becomes 
m=2k) | vas hy — GO] | (x+y?) dx dy = 
PR KR 


The first integral may be evaluated using the equipotential plot and the 
second analytically. Thus the complete state of stress and the torsional 


Moment may be obtained from the electrical analogy. 













CHAPTER II 
MATERIAL AND TESTING TECHNIQUE EMPLOYED 


de Specimen material and dimensions. 

Aluminum alloy AL - 61 was chosen for specimen material because of 
its low shear modulus and ease of machining. A low shear modulus was 
desirable in order to produce relatively large strains and anzles of 
twist of the specimen for a given torque, both of which tended to reduce 
the relative error of the results. A three and one-quarter inch diameter 
shaft was chosen. From Fig. 5 it can be seen that the stress pattern 
must change in going from the small chuck diameter to the relatively lar- 
ger section, this necessitated having a long specimen in order that the 
ends of the ten inch gage length be approximately three diameters from 
the section change. By applying Saint-Venant's principle the stress 
condition in the gaze length could be considered as uniform. 

6. Method of machining specimens. 
The specimens were machined in Seco ndawee with Fig. 5 and Table l. 
TABLE 1 


SPECIMEN DIMENSIONS 





x ~ Ww b 

0° 1.627 in; 332 Sonar, 0 
0 1.576 in, 1. 545 -om. 0.602 in. 
260" lg 57 oats | 1.428 in, 1.331 in, 
46.0° Weg foe bs 1.095 in, 2.265 ime 
65.7) 1.576 in, 0.646 in, 1.436 in, 








a ; = i ff 
| / . 











The specimens were first machined round and then the flats milled. 

All the specimens were checked for axial location and parallelism 
of the flats and were found to be within one thousandth of an inch which 
was considered satisfactory. 

Te Strain gage selection and application. 

The selection of the proper strain gage was an important considera- 
tion as it was to be used in the determinatiom of the maximum stress on 
the flats. The strain varied from zero at the comers to a maximum at 
the center of the flats, this dictated the selection of a gage as small 
as practicable. The principal strain directions are at an angle of 45° 
to the axis of the bar for torsion so a single element gage was chosen. 
Due to the 45 orientation, the width of the gage element was equally as 
important a consideration as the length. Two gages were used, one on 
each flat oriented so as to be ninety degrees to each other, this enabled 
the separation of unwanted bending strains introduced during testing. 
Baldwin SR - 4 strain gages were usede 

A-19 gages with a one-sixteenth inch by one-sixteenth inch gage 
element were initially chosen. The use of these gages was abandoned 
since comsistent results could not be obtained when they were tested on 
the fully round specimen (x=0). Table 2 gives the results obtained 


from an A-19 gage. 
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TABLE 2 


A-19 STRAIN GAGE RESULTS 


Date | Total Drying Time At Total Drying Time Apparent 
Room Temperature At 160-170 F Strains 
4 March 72 hours 0.0123 in/in 
in=-lb. 
6 March 120 hours 040134 1 in/in 
in- e 
12 March 264 hours 20 hours 0.0150 m infin 
in-lb 
18 March 408 hours 62 hours 0.0154 n in/in 
in-lb 


Two separate A-19 gazes were tried but the same phenomenon persisted. 
Fige 6 is a typical curve obtained from an A-19 gaze, the relaxation on 
unloading was prevalent on all runs. In view of this, A-7 gazes were 
substituted and they gave consistent results. The A-7 gazes have a 
one-fourth inch gage length and three thirty-seconds inch width. Two 
AX=-5 gages were installed on the fully round specimen for preliminary 
testing and calibration of the single element zages. All gages were in- 


stalled and dried using the technique recommended by the manufacturer. 
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The maximum angle of twist was less than = of 1 degree so an 











cptical means was chosen to measure angles of twist. The schematio 
diagram for the optical arrangement used at each end of the 10 inch 


gage length is shown in Fig. 7. 


collimating lens 


tp hair line 3alvanometer 


VS lens 
filament image 
filament | 
O 
@ | 
mirror 
hairline image 
wax 


frosted scale 


Fize 7- Schematic diagram for optical system used for angle of 

twist moasurement. 

The essential elements are a hairline illuminated by a light source 
so arranged that the rays from the hairline are reflected by a mirror 
mounted on the specimen and focused by a galvanometer lens, on a frosted 
scale. The actual equipment used can be seen in Fig. 8. The hairline 
images from both the light sources were focused on one frosted scale. 
Thus differences in the positicns of the images on the frosted scale was 
a measure of the anzgle of twist in the gage length. Since small angles 
were involved, the sine, tanzent and anzle in radians were considered 
equal. The formula for the angle of twist for a given load is, 


p = oo where L\h is the change in distance between the hairline 


images and d is the mirror-to-scale distance. The angle of twist per 
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unit length for a given load is e= 7% where 10 inches is the gage 


length. 
9- Strain gage tests. 

The equipment as used is shown in Fig. 8. The specimens were 
loaded in a standard Tinius - Olson Torsion Tester on the 0 to 30,000 
in-lb range except for the thin specimen ( Cc = 65.7) when the 0 to 
12,000 in-1lb range was used. Two AX - 5 strain gages mounted on the 
fully round specimen at physical locations ninety degrees apart were used 
to determine the best method of putting specimens in the machine without 
introducing extraneous loadings. If the rotating head of the machine 
was at or near zero degrees and the chucks were tightened slowly and 
alternately the machine applied essentially pure torsion. Two runs were 
made on each specimen, one with each gage in the up position. A Baldwin 
Southwark Type K Strain Indicator was used with a Baldwin S R- 4 12 
Channel Switching and Balancing Unit to measure all strains. fhe atl! 
erature compensating gages were moumted on a separate piece of the same 
material. The specimens were loaded and wmloaded in equal increments. 
After each increment of load had been applied the strains, hairline 
image positions om the frosted scale and load were recorded. Strains 
were read to the nearest [4 infin. load to the nearest 10 inch-lb. 
end hairline image positio to the nearest tenth of a millimeter. 

10. Analog Experiment 

The equipment used for the analog experiment is show in Fige 4. 

It is a General Electric Analog Field Plotter. The experiment was car= 
ried out mly for x=45, woe3O /n. As can be seen in Fig. 3, 
Cx is the angle between the x axis and the exterior comer of the 
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Fige Ge Experimental Equipment Set Up 
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of the cross section the conducting sheet represented only 4 of the cross 





section as shown by the shaded area in Fiz. 3. The large dimensions 













were chosen to increase the accuracy of the results. 

The null method of determining potentials made it convenient to 
record experimental potentials P on an arbitrary scale. Accordingly, 
the potential at the center of the flat (w ,0) was taken as O and that 
on the curved portion of the boundary was taken as 1. As is shown late 
the potential V of Eq. § is related to the experimental potential P by 
the equation V=1+P. 

The parabolic potential distribution along the boundary X = W was 
approximated by a series of linear potential dropse Fig. 9 shows the 
desired parabolic distribution and how it was approximated by a series 
of linear drops. There was no net gain or loss of area wmder the curve. 
Small increments were taken at the center of the flat in order to have 
the smallest deviation from the desired distribution, and keep errors in 
potential gradient to insimnificant magnitudes in this region. Table 3 


gives the potentials maintained at the chosen locations along the boundary. 
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TABLE 3 








POTENTIALS MAINTAINED ALONG BOUNDARY 


X (ins) P 

0 0 

® 0 

2 0.0039 
4 0.0162 
7 0.0496 
12 0.1554 
17 0.3165 
22 0.5332 
| 0.8054 
30 1.000 


The experimental potentials measured along the boundary are plotted 
in Fig. 9. To establish the above potential points along the boundary with 
linear drops between points, 0.004 inch diameter Advance wire with lead 
wires soldered to it at the selected locations was clamped onto the 
boundary of the paper by the clamping arrangement shown in Fiz. 4. The 
lead wires were cammmnected to a slide wire potentiometer which was across 
the output voltage of the rectifier wit. The actual bowmdary potentials 
were established 0.15 inch in from the wire by varying the positions of 
the lead wires on the slide potentiometer. It was found that if the 
correct potential distribution was set alony the boundary wire the act- 
val potential distribution in the adjacent paper was significantly dif- 


ferent. The comstant potential alons the curved portion of the boundary 
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equipotentials were plotted using standard null technique. 


















CHAPTER III 


REDUCTI CN OF DATA 












1]. Strain sage calibration. 

The AX-5 strain gages on the fully round specimen were used to 
calibrate the single element A-7 gages. This was necessary beoause the 
gaze factor given by the manufacturer gives correct results for uniaxial 
stress in the direction of the gage length. The gage factor for this 
strain condition was fowd from 


ae A-7 Manufacturer's gaze factor 
~ | A=-7 gage factor for this strain condition = 


where € is the strain in the direction of the A-7 gage as obtained 
from the AX-5 gages and © is the apparent strain obtained by the A~-7 
sazse when the manufacturer's gave factor was set on the indicator. 

The manufacturer's gage factor was 1.95 and the apparent ga:e factor for 
this strain condition was found to be 1.97. It is noted that this re- — 


sult implies a negative transverse sensitivity for the A-7 gase. This 


is in agreement with transverse sensitivities reported by Baumberger 


and Hines Lae. 
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e center of the flat and 





the computations for torsional moment were carried out as discussed in 
section 4. By Eqs. (8) and (OS) , P=K(V) = ee (27) 

At (wu ,O) j KV, me ae) and at (ur, w) , K\%=G Gea 
Takinz the applied potential difference as wity, 4Y-Y4e=l 


E» Ogee ” = 
Then K= 3 : y= 2m , aa 
point on the surface is, V=V+P=/+P where P is the expe 


- The potential at any 


mental potential. Fig. 10 shows the equipotentials over a quarter of 


the cross section. At (ur, O) Tyg =O and 
T= 8% +G0x =-SBe9P) it Gew 
id ox a XJ (ur, O) 
(olisak is obtained from Fig. 11. Eq. (4a) now becomes 
0X} (ur, O) 


M, = @@-ur* [Jf I+P)A&X dy — oo] Jotun 


Simpson's Bega rule for approximate eae was used to obtain 
the volume wder the potential surface for evaluating the first integral. 


The value of the last integral was obtained analytically. 
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Fize 10. Equipotential Traces of blectrical Analog. 
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RESULTS AND DISCUSSION 


13. The results. 






Numerical results obtained from tests are given in Table 4. The 
















values for x=45 were obtained by the electrical analogy. 


TABLE 4 
VALUES OF MAXIMUM SHEAR STRESS 


AND TWISTING MOMENT 


xX Bi Ml, 

0° 1.00 6Q@-ur 1.56 GOr?* 

ii 1.13 GOu 1.5760rt 

25° 1.31 GO w 1.44 GOr* 

45° 1.66 GOW 0.964 GOr* 

46 1.56 GO Uu 0.923 GOrt 
’ 65.7" 1.61 GG uw 0.257 GOr* 


rm is the radius of the shaft. 
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14. Comparison with 


Fig. 12 gives a comparison of the present results with those of 


previous investigators. Carter and Oliphint stated that they believed 


their results were accurate to within 2 or 3 per cent. Okubo gave no 


figures for accuracy of his results. 
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15. Sources of errore 


The strain vs torque relationship obtained with the AX-5 strain 
gages on the fully rowd specimen was linear within reading accuracy. 
The extraneous bending stresses introduced by the torsion machine were 
found to be less than 1% of the torsional stress. Results were cor- 
rected for bending stresses by averaging readings of gages on opposite 
flats. The orientation of the strain gages is believed to have been 
within l°of the desired 45° orientation. Location of the mirrors on the 
10 inch gage length was within + 1/16 inch. The strain gage elements 
covered a finite area on the flats of the shaft and the stress varies 
from zero at the exterior corners to a maximum at the center of the 
flat. Thus the strains indicated by the gages were not the maximum but 
some sort of average over the gaze element. This error is believed to 
be small. 

The manufacturer of the Teledeltos paper states that it has aniso- 
tropic variations of tz to £8% maximum deviations from a mean value of 
resistance. The rated linearity of the potentiometer used to establish 
the equipotentials was the No correction was made for the boundary 
potential being established 0.15 inch in the interior. 

The optical system used for measuring angles of twist gave reading 
accuracies of 1% of the total angle of twist. The obtaining of greater 
acouracy would have greatly complicated the optical system. 


The overall results are believed to be accurate to within 4%. 
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CON CLUSI ON 
The experimental results checked closely with the previous 
mathematical solutions to the problem. Althouzh the results obtained 
by the electrical analozsy were reliable, the work involved was many 
times that of the direct experimentation for the same specimen 


configuration. 
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